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THE DIFFERENTIAL EQUATIONS OF DYNAMICS 

BY 

AETHUE C. LUNN 

§ 1. Introduction. 

The general problem of the dynamics of a system depending on a finite 
namber of parameters may be considered equivalent to the problem offered by 
the canonical differential equations : 

dx, .dF dy, dF 

indaded in the general type of equations in the normal form : 

(2) ^ = X,(a;....«„0. 

With respect to the properties of the solutions defined by equations of this type 
most investigations have naturally referred to the case where the functions X. ^ 
and consequently the solutions x.^ are holomorphic with respect to all argu- 
ments, and even where real variables only are primarily considered, admit direct 
interpretation for complex values. The researches of Poincab^ in astronom- 
ical dynamics,* which have contributed a body of theorems of great importance, 
have been based mainly on this assumption, though many of his results admit 
an obvious extension to other cases. 

A namber of writers have however found it possible, by restricting the varia- 
bles to real values, to secure a variety of results with more general conditions 
on the functions X^^ . In particular, dynamical applications of the method of 
real variables are either developed or directly suggested, for example, by papers 
of PoiNCAR]&,f Bekdixson,^ Picard,§ Painlev^, || and Hadamakd.^ In 

^Colleoted for the most part in Le$ MHhode» NcuveOea de la MSeamque C^Utie, Paris, 1892, 
1893, 1899. 

t /9ur lea Oourhea diflnies par lea Squationa diffhewLidUa: Journal de Math. , 1881, 1882. 

t Same title, Acta Math., 24, pp. 1-88. 

§ 8ur ViniSgration de eertainea S^ii4ma dijfh-etUidlea par lea aSriea: G. R, 124, pp. 214-217. 

\\8ur lea mauvementa et lea trajedoirea riela dta affaUmea: Bull, de la Soo. Math, de 
France, 22, pp. 136-184. 

^8ur eertainea prcpriiUa dta trojeeioirea en dpuunique: Jonr. de Math., s^r. V, 3, pp. 

331-387. 
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4 A. C. LUNN : THE DIFFERENTIAL EQUATIONS OF DYNAMICS 

the present paper certain known theorems relative to equations of type (2) are 
somewhat extended and the results applied to problems of dynamics, as repre- 
sented by systems of type (1) or by certain auxiliary systems, attention being 
confined chiefly to the case of real yariables, as arguments of functions not in 
general supposed to be analytic. 

The case of chief interest for dynamics, perhaps, is that where the solutions 
x^ = <f>^{t) are uniquely determined by the initial conditions <^j(^o)~^{o) 
though sometimes other determining conditions may be important, for instance 
that a system move from one assigned configuration to another. These solu- 
tions, as functions of the initial values x^^ t^^ and of parameters which may 
appear in the differential equations, have been studied, with special reference to 
continuity and differentiability, by several writers, in papers which are referred 
to later.* Some of these results are here deduced as corollaries from a more 
general theorem on the behavior of the solutions as functions of a parameter 
taking the values in an arbitrary point-set with a limit-point, a theorem also 
shown to have as corollary an extension to differential equations of the term-by- 
term integration of infinite series. This is illustrated by the analysis of reso- 
nance in a vibrating system, and by certain methods of approximation in celestial 
mechanics. In particular there is given a set of conditions sufficient to ensure 
that the successive approximate solutions furnished by the method of '^ variation 
of parameters" shall converge uniformly toward the exact solution during a cer- 
tain interval of time, together with a proof that in a very general class of cases 
this method can be defined so that these conditions are satisfied. 

§ 2. The Cauchy-Lipschitz Approximations. 

As basis for the following study of the solutions defined by the differential 
equations it will be convenient to use the convergent approximations employed 
in the existence-proof of CAUCHYf and Lipsohitz,^ partly because of the 
important property to be stated shortly. For the solution x, (^), • • *, a;^(^) or 
P(f ), starting, for f = ^q, at the point P^(ajj^,, • • •» a^o)' ^^^se approximations 
are obtained by partitioning the interval t^- "thj numbers in order f^, t^^ •• *, 
t^=t^ and successively defining 

Because of the dependence on the set of numbers marking the partition p , we 
shall denote such an approximation by a;j(/>, ^), or P(p^ t). Then beside 
serving to prove in a general case the existence of the solution P(f ) for a cer- 

* Ree Dotes to § 5. 

t MoiGKO : Legons de Calcul. • • •, II, p. 385. 

} Lehrbuch der Analffsia, II, p. 600. 

§ G. R, 126, pp. 136a-1366. 
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tain ^interval, these approximations, as shown in papers of Pigard§ and Pain- 
LEV^,* possess the important property stated in the following lemmas : 

Let Z> be a closed realm of points P{x^j -"t x^)^ and D the complementary 
realm of points JP{x^y • • •, S^), and t^-- - T b, closed interval, where T may be 
greater or less than t^. Let X.(^x^^ ••*, x^, t)^ i^l ••• n, be defined as 
single- valued functions, for (a;,, •••, x^) any point in D and t any in the 
interval t^-- - T^ with the following properties, provided (a|, . • •, x'J, (x[\ • • •, ac^') 
are in -D, and ^, t\ f in ^^ • • • 2^: 

A. There are n positive numbers m^, • • *, m^, such that 

\^i{^[.-'.<.t)^X,{x';,.^-,x:,t)\^m'f,\x]-x';\. 

B, Any positive number e has a corresponding pair of positive numbers 
5(€'),S(€'),8uchthatif |^'-f"|<^(€'),and|»;-a;;/|<S(€'),j = l ..-n, 
then 

(5) I X.(«;, . . ., x'^t t') - X,{x\\ . . ., «;, n I < e'. 

Let M be the finite upper bound in Z> and t^- --Toi the moduli | ^^ | , which 
exists because of (B) : let e, €, be two positive numbers, and ^^ • • • ^^ + r a closed 
interval, where r lies between zero and T— t^i define 

(6) A* = Z^i. ^^°2^ e^^-l) ^ 7(t, c) = smaller of 5(6), -gr-. 

Lemma 1. ij^ «,(<), • •-, x^(t)^ or P(^), wJiere J^it^) = Po, exists as a 
solution of the differential equations (2), over the interval t^-- 't^+ Tj and such 
that for all points of D and all values of t in that interval the lower hound 
of the sum 



is not less than e + ei then for a partition p toith every interval less than 
7(t.€): 

(7) T,\x,it)-x,(p,t)\^, 



and for any two such partitions p\ p"\ 






(8) 5:i«,(y>o-a'.(p".oi<«- 



Conversely, with r as above, if P( j? > ^q + t) exists as such an approximation, 
with intervals less than /(t, e), and is such that for all the « + 1 points 

*C. R., 128/ pp. 1505-1508) and Ball, de la Soo. Math, de France, 27, pp. 149-152. 
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(«,., • • • , a5„,),^' =s • • • «, oocurring in equations (8), and for every point in 2), 
the lower bound of the sum 

is not less than e + c, then for t any in the interval t^ -" t^ + r let p^ p\ p be 
partitions with intervals less than /(t, 6). 

Lemma 2. Under these conditions the partitions p , p\ p' define points in 
the realm -D, and the approximations P{p\ t)^P(^p\ t) satisfy the inequality 
(8). There is a solution P(^t) of the differential equations^ lying within the 
realm over the interval t^-'-t^ + T^ satisfying the inequality (7), and such that 
P(fj,) = P^. For all points of this solution in that interval and all points 
of D the lower hound of the sum 



n 



1 

is not less than e . 

In particular the conditions of these lemmas are always satisfied if t, as above, 
is besides not greater than g/Mj and ^ + e is the lower bound for all points of 
D of the sum 

n 

1 

The Cauchy-IApschitz approximations con/oerge^ and the solution is unique^ 
as long as the points P{t) are within the realm D.^ It is here understood 
however that the computation can be thought of as begun with intervals of any 
required fineness. 

§ 8. The solutions as functions of a parameter. 

With the aid of these lemmas the solutions will now be considered as func- 
tions of a parameter r, which takes all values in a point-set R admitting r^ as 
limit-point. The set R is to contain all its limit-points except r^, so that the 
addition of r^ to the set gives a closed set which may be called R^. We shall 
use the term ^^ neighborhood of r^ " to denote all points of R lying within a cer- 
tain distance of r^, and denote it by jS( ), with the quantities on which that 
distance depends indicated within the parenthesis. If jS is not bounded infinity 
may be thought of as one or more definite limit-points, according to convenience 
in making special interpretations of the results. 

Over the realm D and closed interval 2^^ • • • 2^, where T^ < T^ let 
X\(^x^^ * • *> ^n' ^)' ^^^ ^^y value r in jB, and X^(^x^^ •• •, jc^, f), associated 
with the limit-point r^, be defined as single-valued functions with the following 
properties : 

*Painlev£, PicABD, loo. oit.; lee also LindbIi5f: Jour, de Math., s^r. 5, vol. 6, pp. 
423-441. 
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I. \X^,{x\t)-X:(x\t)\mm,Y.\x;^x';\ 

where m^ is independent of r. 

II. The functions X"^ are continuous in D and T^-^ T. 

III. For any positive number e there is a neighborhood ^(c) of r^, such that 

\X,{x,t)^X\{x,t)\<€ 

if r is in R{e) • If R has a limit-point /, distinct from r^, it will also be sup- 
posed that the corresponding function satisfies a condition analogous to IH, with 
X^ replaced by X'^ . 

The functions X. are therefore continuous, being uniform limits of continuous 
functions. Moreover, since the realm 2), interval T^*" 2^ and set R^ are closed, 
a modification of the familiar theorem of uniform continuity yields the following 

Lemma 8. For any positive e there are the positive numbere S(€), d(e) 
depending only on €, and such that if \x'j — a5J'| <S(e) and |^' — f\<0{e) 
then 

(10) \X{x\l!)-X{x",f)\<€ 

where X is any one of the Junctions XT, X.. 

For the contrary would imply the existence of an infinite sequence, determin- 
'ing a function, either X^ or an X*[, and a point (aSj, • • • , sc,, t), such that in 
any neighborhood of that point would be two points for which the values of that 
function differ by at least € , a violation of continuity. The moduli \X^\j \ X\ \ , 
have thus a single upper bound M. 

Now if (a5(, • • •, aj^), (a|', • • •, aj^'), be two distinct points in 2), t a value in 
7^ ... 2^, and n^ any positive number, there is by III, I, a function X^ satisfy- 
ing the following inequalities : 

whence by addition : 

\X,{^,t)^X,{x\t)\<{m,+ 2n,)Y,\x]^x]'\. 

i 

But n^ is any positive number, so that this is equivalent to 
Lemma 4. The fwivetions X^ also satisfy the condition : 

(11) \X,ix\t)-X,ix'\t)\^m,J^\x'j-x^\. 

J 

It may be noted for example that condition I wiU be satisfied if the functions 



8 A. G. LUNN : THE DIFFERENTIAL EQUATIONS OF DTN AIHCS 

X], X^ have continuous partial derivatives dX^jdXj^ dXJdXj satisfying rela- 
tions analogous to III above. 

We now compare the solutions sc^(^) or P(^) of the equations 

dx 
(12') ^*-X,(a!.,...,a!..0 

for the initial conditions P(t^) = P^(a3j^, • • •, a;^^,), with the solutions x'l(t) or 
P''(^) of the equations depending on the parameter r : 

dx^ 
(12") -^-Xj(«r,...,a;:,0 

for the initial conditions P'"(^J) = Pj(a5jjj, • • •, fic^o) » ^^^^^ ^^ initial values 
»^^, tl also depend on r in such manner that 

IV. For any positive numbers e\ e\ ff^ there is a neighborhood R{e\ e\ ff) 
of r^, such that if r is in the set R{e\ €\ ff)^ then 

(13) Zla^Jo-^^K^'' \ti-h\<ff. 

Now let t^ be interior to the interval 2^ • • • T, r a positive number not greater 
than T^t^^ and e, €, any two positive numbers. 

Theorem I : If P(t) is a solution of {12') for all values of t over the interval 
^o" '% + ''' 9 ^^^ ^^^^ ^^^^ /^ ^^ those values of t^ and all points of D, the 
lower bound of the sum SJ ^^ — ^i(^) I ^ ^^^ less than e + e; then there 
exists a neighborhood Ii{e) ofr^^ such that if r is in the set Ii{^)^ equations 
(12") have a unique solution P'^t)^ determined by the condition P'"( ^J) = jPJ, 
defined over the larger of the intervals tQ-^t^ + r and ^J • • • t^ + t, and such 
that for every value oft in the common interval 

5:|a5,(0-xI(OI<e. 

The lower bound of the sum 2^ | x^ ^^\{t)\is not less than e. 

Proof Let the equations (8) be formed with the same partition p for equa- 
tions (12^) and (12"). To fix the ideas suppose ^^ < t^, and the interval ^^ • • • ^o 
partitioned into k intervals by t_j^ = ^J, ^-jh-p • • • » *-i » ^o* Then equations (8) 
may be written 

^ij+l — ^ij = ^i{^<rJJ tj){tj^i--tj) (J = 0...5— 1), 

^i,j^h+i — a?,i+* = ^i{Kj+k9 h){^j+t — tj) (i=—*, • •, 8—1), 
whence 

(«*j+i--a^ij+*+i) = (a^u--<i+*)+ {-y<(«^.i, «i)--^i(»..i» ^i)} 

+ {x:{x,j,tj)^x:{x,j^,,tj)}{tj^,^tj}. 

With the abbreviation 
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(14) Ei««.>-"?.i+*i =>»;(;'). 

i 

we have therefore by (9), I, and (11) : 



i 



With the positive nnmber X to be defined pat 

(16) sKp)+l=yKp)^ 

then by (9), m, if r is in the set ^(X/m): 

'^;+i(i')-^<[^;(;')-^] + {'^+/*[^;(;')-^]}{<i+i-<,} 

whence 

T7+i(i')<^;(i'){i +/*(«*+!- <i) }< ^;(p)«*<^*-^' 

and by multiplication 

(16) F;(|>)< FJ(;>)«''(^*" (i=l.--,.). 

Bat 

F-(2>)-^ + El«...-«i.*l 

where 

2:|a!«,,-a«.*l = Zl»'i.o-a?.oH-ZI«^*-'«^«l 

i < < 

also 

Zi«^*-*f.«l=»-^-i«o-«;i 

BO that (15), (16) give 

'S'5 (i>) + ^ < { ^ + »^- 1 «, - <J I + E I «*. , -« J. , 1 1 «"<'>-'•' 

or finally by (14) 

issn \ 

(17) s l^X/*. «)-a?(i'. Ol<{»-a''K-«;i+Ela'«-'«s.l}«'''*+ ii {«"'"- 1 

where ^ is any in the interval ^q ** • t^ + r, and r^ » T— ^. The same form 
would result for the interval %-• -t^ + t from the supposition tl>t^. 
Put now 

(18) 

^as smaller of g—j^e-"^', <, — T; 

then if r is in the set Ii'(e', t\ ff) defined by the context of (18) : 
(19) E|a'i(l>,0-«<(i'.OI<^- 
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Moreover this single partition p, by (10), (11), and Lemma 1, can be taken 
with intervals less than /(r^, e/6 ) defined by (6); in which case 

(20) El«i{«) -«*(/>. 01= I 

and by combination of the last two inequalities : 

(21) 2:ia'.(0-»:(p.0i<|-- 

Inspection of the preceding inequalities shows also that the solation P(^) or the 
approximation P^{p^ t)y according as ^^ > or < ^jj, can be carried back so as 
to refer to any point in the interval <^, — ^ • • • ^^ + t, the inequalities (19), (20), 
(21) remaining fulfilled ; and the lower bound of the sum ^.\x^^x^(t)\ would 
not be less than e + |«- Also the lower bound of the sum 2^ | ^^ — a;^ (/> , t)\ is 
not less than e + ^e. By Lemma 2, with € replaced by c/6, the conditions for 
the existence of the solution P^(t) are fulfilled for the interval ^^ — ^ • • • ^g + t, 
and the theorem is proved, since the combination of inequalities (20), (21), with 
the inequality 

(22) 2:i<(o-<(/'.oi^g 

gives 

(28) El^iCO -«;(«) l<e. 

t 

By a process so closely parallel to the above that it need not be given here 
can be proved the analogous theorem : 

Theorem II. ^ there is a point r^ in the set It{e\ e', ff)y defined accord- 
ing to (18), and for some point t^+ t in the interval t^"- Ta corresponding 
approximation P'ipy % + t), with partition-intervals less than I{T^y ^/6)» 
for every point of which and for every point of D the lower bound of the sum 
^^\x^ — a^^l is not less than e + e; then ifr is any point in R{e\ e', ff) and 
p any partition with intervals less I^r^^ ^Z^)* ^ solutions P{t) and P^(t) 
and the approximations P{py t)^ ^{p^ ^** <w^^ ^he interval t^^ff "'t^+'^j 
the inequalities (19) • • • (28) are satisfied ; and for any of these solutions or 
approximations the lower bound ofS^-lx. — ficj is not less than e. 

In a similar way the solutions corresponding to two different values of the 
parameter can be compared with each other instead of with the limiting solution 

The above theorems may be considered a direct extension of the existence 
theorem of Caught and Lipschitz, in that the existence of a solution for one 
system of equations is used to show the existence of a solution for another 
'^ neighboring " system. Several special forms are at once suggested. 
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If r is a point (/li^ • • •, /x^) of the ^^^mensional continuum in the neighbor- 
hood of the point r^, the conditions given are sufficient for the continuity of the 
solutions as functions of parameters in the ordinary sense, in particular as 
functions of the initial values x^y t^. This case has been treated repeatedly 
by previous writers, together with the question of the existence of derivatives 
with respect to the parameters. This latter question will be taken up below in 
§ 5 as a direct application of the preceding theorems. Of more direct interest 
for the sequel, however, is the application given in the next section to differential 
equations whose second members are infinite series. 

§ 4. Equations with second members infinite series. 

The following iUustration of the results of previous sections gives a direct 
extension to differential equations of the theorem that a uniformly convergent 
series of continuous functions can be integrated term-by-term. 

Let r^ be the definite point at infinity, and let the set li include all points 
(iVj, • • •, iV]^), where the IPs are positive integers or zero. Then the X^'s 
may be thought of as defined by the finite series 

(24) Xj = 2Xf» "*(«!, •••,«,, (m«=O...JVa,a=l...ifc) 

and the X/s by the infinite series 

(26) X^=r2X*> "**(«!» •••5 «H> («la = 0- •• 00, a=l •••*). 

Then the conditions imposed in the preceding section are sufficient in order that 
the solutions corresponding to the finite series shall converge toward those cor- 
responding to the infinite series. An example of considerable generality is given 
in the following theorem, where for simplicity we suppose x\^ ^x^.t^^t^. 

Theobem III. If the functions -2r7^""'"**(a5j, •••, ac^, t) have partial 
derivatives with respect to a;^, • • '^ 90^$ which are continuous over the realm D 
and interval T^^^T^ and the infinite series (25) and 

d 

(26) Z;^.= 2^Xf ** (iii« = 0.-.oo,a==l...ife), 

converge uniformly there; if moreover the equations whose second members are 
defined by (25) have a solution P{t) lying iiUerior to the realm D at every 
point of the closed interval ^q • • • ^q + t, then there is a point {N[ » • • •» -^i) 
such that for every point (iVj, • • •, -^), where N^ = N'^^ the solutions of the 
equations 

cfojT 

(27) -^-SXf -'(ajj, •••,<,«) (in.=0,.--,iV., a=l,..,fc), 

exist and lie interior to the realm D oner the interval t^---t^-\- t; and at the 
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point (iVj, •••, iVj) approaches infinity those solutions converge uniformly 
toward the solution P{t), 

For the functions JT^ are then continuous, and have continuous partial deriva- 
tives, which are the functions X^ of (27) ; * the moduli | X^ | , ^' =» 1 . . . n, for 
the finite and the infinite series have an upper bound m^, which may be taken 
as the factor so denoted in (9) I ; and since P(^) is interior to the realm D 
over a closed interval the lower bound of the sum SJ^^ — aSj(^)| is positive. 
The conditions of Theorem I are therefore satisfied. It is obvious that the state- 
ment of Theorem II, also, presupposing the existence of some jP'(^) instead of 
P(^), can easily be adapted to this case. 

If the interval ^^ • • • ^^ + t is sufficiently short, for example less than gjM^ 
where g is the lower bound, for all points of ^, of 2J a?^ — x.^ \ , and M is the 
upper bound of the moduli | XJ , the solutions all lie within D over that interval, 
and the point (iV|, • • •, N*j^) may be taken as (0, • • •, 0). 

A simple case of interest is that of the linear equations : 

(28) ^* = i;^<,(Oy, + V^<(0, 

for which the general solution has the form : 

(29) y.-EVOJy^ + EjT'^rf*}. 

the functions 4>^^'^^j being continuous; in which 0^(^t^)^ 8^, and A is the 

,(i,^'ssl...n), with minors A^ defined by the conditions 



determinant 1 0^ 






where we put 8^^ = or 1 , as a 4^ /8 or a = /9. Then the preceding theorem 
shows that if the functions 0, y^ are uniformly convergent series of continuous 
functions, the functions ^, their derivatives 6\ and the integrals in (29), can be 
written as series of the same character. 

An illustration of this is found in the theory of resonance in vibratory 
systems, which leads to equations of the type (28), in which the <^'s are con- 
stants, and the '^'s represent the external forces producing the ^^ forced " vibra- 
tions. The theorem then gives conditions which justify breaking the particular 
solutions into series, each term of which corresponds to a component of the 
disturbing forces. If the external forces can be represented as uniformly 
convergent series of components which are continuous functions of the time, 
then each component produces its own distinct effect in the resonance of the 
system. 



See 6. g.: Jordan, Cours d^ Analyse^ I, §§ 328, 330. 
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It is obvioas that in the general case of Theorem III the solutions x^{t) can 
be thrown into the form of uniformly converging series 

(80) a5,=r2/7l **(0 (ffla = 0-OO,a = l...ife) 

where the terms are defined uniquely by the conditions 



m^=sN^ 



(81) 52 fT **(0=»f*' ••^•(O- 



»a*0 



Combined with Lemma 1, this theorem would yield, for example, the representa- 
tion of the solutions by infinite series of polynomials in ^, as deduced from the 
Cauchy-Ldpschitz series by Pioard * and PAiNLEy:ig,f who do not however make 
clear the fact that the multipliers m. in Condition I above, or its equivalent, can 
be considered the same at every step of the process, in other words independent 
of the parameter r, a property apparently vital to the proof of convergence. 
That this could not be proved from II, III, and the milder form of I where m^ is 
replaced by m^ , appears from the fact that not every continuous function satis- 
fies Condition I, but every such function is the sum of a uniformly convergent 
series of polynomials, each of which satisfies that condition. In the case in 
point, however. Condition I in the form given above is satisfied, provided the 
polynomial approximations are so chosen that they converge uniformly to the 
^.'s, and their first partial derivatives converge uniformly to the partial deriva- 
tives dXJdXj. 

§ 6. The avxUiary linear equations. 

With the notation of § 8, let ^(r) be a single-valued function, defined at 
every point of the set R^ and to avoid trivialities let the modulus |0(r)| have 
a lower bound distinct from zero. Then equations (12'), (12"), give : 

(82) $-i^i(f:.--.f;.o 

where 

(88) fj-<^(r){x:(0-»,(<)}, 

and the functions Y\ are defined by 

(84) rj(y.,...,y„,0=-^(r){x:(»+|,<)-X,(a!,«)}. 

It is understood that x^'-'X^ are thought of as the functions of t defining the 
solution P( ^) of the equations (12^). We now add to the hypotheses of Theorem 
I the following conditions : 

*C. R., 128, pp. 136^1366. 
tC. R., 128, pp. 1605-1508. 
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V. For every set of positive numbers e", ^\ ff\ there is a neighborhood 
R\e\ e'\ ff') of r^, such that if r is in that subset, then 

(86) \^{r){X',{x,t)-X,{x,t)]-X\{x,t)\<^' 

and 

(86) |*(r){«:,-xj-f:.|<e", |*(r)(«;- «.)- '».l<<^ 

where A^ (as, ^) is a continuous function of a?^ • • • x^, t over D and 2^ • • • T, and 
fio» '',i'no^%9iQre8X numbers. 

VI. The functions X^ have partial derivatives 

continuous over D and T^-*- T. 
We may therefore define 

(87) F,(yi, • • •, y., = -^i(«i» • • •» ««» + E-^</(^i' ' ' '» ^..^ OVr 

Theobem IV. Under these conditions there is a neighborhood R'^ of r^, 
suoh that ifrisinR^y then 

ZlfJ(0-f,(«)l<« 

where the ^^{t) are defined by the equations 

(88) §-^;(«.o+E^<^(*.of, 

with the initial conditions 

The proof is immediate, since the functions Y\^ Y^^ satisfy conditions precisely 
parallel to I, II, III, above. 
Condition I: 

rj(r,o-^:(r,o={^*(«'+f.«)-^j(«'+|.«)}-* 

so that 

j':(r,o-J":(f.«)i^»»,i*iE|^^'|-'»i2:if;-f-: 

ik I T* I It 

Condition 11 : Y\\% continuous since X\ is continuous, and ^(r) is nowhere 
zero. 

Condition III: the difference Y\{i^ t) — l^({i "^7 ^ written as the 
sum of the expressions 

4>-{X\{x + il>i>,t)-X,{x + il4>,t)}-X\{x + il^,t) 

and 

^ • { X,(x + f/<^, t) - X^{x, t) } - 5:^<,(a!, Ofy 

i 
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and 

X',(x + tf<l>,t)-X',(x,t). 

The first of these is controlled by V, the second by the theorem of finite differ- 
ences,* the third by the continuity of X\ . That the I's may be considered as 
lying in a bounded reahn appears if both sides of the inequality (17) are multi- 
plied by 1^1, and the right-hand sides inspected with the conditions imposed in 
mind. As to the initial conditions we have : 

f:(«.) = *-{a':(«o)-'«.(<,)}-<^-{<(«»)-«:(<s)} +*•{<(<;)-«'*(<,)} 

or 

^i(%) = <!>■ («« - «'«) -<t>{ti- <o)K(<.) - «:(«;)] /(«. - «;) 

which may be compared with (86) and (89). 

If the set H is the linear continuum, as realm of the single real parameter r, 

and ^(»') = l/(»' — ^o)» ^^^ f^ *^® ordinary derivatives, and we recover the 
results obtained in the memoirs cited in the footnote.f The derivatives dx./dr 
are defined by the equations 

(40) . ^/^_?^..2X.^ 

with the initial conditions 

For example, if we put X^(x^^ - - *, a?^, t) « <f>^{t)y then the equations 

dx. dx. 

define n + 1 solutions of the homogeneous linear equations 

(42) |'-i:*<,(Oyr 

There must therefore be a linear relation between them. The initial conditions 
are 

Mfe),.,.- k - '•• (I'),...- - ^'<^' • • ■• "• ^>- 

* Jobdan, Ooun, I, { 91. 

t For oontinaity aod differentiability with respect to initial values and parameters see 

Pbano: Torino Atti, 33, pp. 9-18. 

NioooLBTTi: Rom. Ac. d. Linoei (5), IV, part 2, pp. 316-324. 

EscHSBiCH : Wien. Ber., 108, pp. $21-676. 

LiNDBLdF: J. de Math. (5), 6, pp. 423-441. 

SoHUB: Math. Ann., 41, 569. 
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We find therefore, for t ss t^ and consequently for every value of f , 

(44) a^' + 2:^*(«'..»---.a'-,.O^-0.* 

With respect to equation of type (42) it may be noted that if y^, X = 1 • • • £, 
be k solutions and we put 

(*5) f<i...i*=-|yij (x,A*=i-.*), 

then 

a system of homogeneous linear equations satisfied by the determinants (45). 
From the case £ s n comes the formula of Liouville : 

(47) D{t) - 2)(e,)e-^^^?*'*^* 

where D is the determinant | y^ 1 9 (t ^J » 1 , • • '9 n). For the system of solutions 
dxJdXjQ we have i>(^^) » 1 , by (43). This system is therefore a fundamental 
system. 

This shows for instance that if the solutions of (12') be written 

then the equations 

(48) x^ + Aaj, -/.(e, e,; x,, + ^x,,, • • •, a;,o + ^»o) 

can be solved to express Aa;jo9 •••9 Aoi;^ uniquely in terms of Ax^^ •••9 Aas^, 
provided the moduli of the latter are properly restricted, since the determinant 
D is the corresponding Jacobian of the functions /, which have moreover 
continuous partial derivatives with respect to all arguments.f 

The same fundamental theorem on implicit functions in real variables shows 
that the proofs given by Poingar^ ^ of the existence of families of periodic 
solutions depending on a parameter apply in part to more general cases than 
those involving power series. 

Suppose the functions ^^(^9 fi; x^^ ••*, x^) have partial derivatives with 
respect to a;^, • • *, as^, fi, which are continuous functions of all arguments, and 
that the solutions of the equations 

dx 

(49) ^'=iE;(«,M; a5„ •••,«.). 



*Compare Bkndixson, Ball, de la 800. Math, de Fr., 24, pp. 220-226; Livdbl5f, 
loo. dt.; HoLMOS£N, Stook. Akad. Bih., 26, part 4, p. 19. 
t Jordan, Coun, I, §92. 
t Mithodes NouvdUi, I, pp. 89-89. 
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supposed written in the form 

(50) «.— /(^Z*; *o»«io^ •••»a'«o) 

are such that for a particular set of Talues ^o> A^o* ^lo* ' "y ^no ^® solution has 
the property 

(51) «««/,(«o+2^'A*o; ^o»«'io»--»»«o) 

then the equations 

(52) /,(^o + 2t,A*o + V; %+At,,x,,+Ax,,,''',x^+AxJ--(x^+Ax^)=0, 

may perhaps be solvable for the increments of the initial values in terms of Afi. 
Two such cases are obvious from the work of Poincab^, without further restric- 
tions on the functions X.. 

First, if the X/s involve t explicitly, with period 27r, there will exist a family 
of constant period 27r for all values of fA in the neighborhood of fbg, if the 
determinant 



(58) 



JO 



{»,i = l--n), 



is not zero for t^t^ + 2'7r. In this case the condition A^^ » is imposed 
because of the explicit appearance of t^ and Ax^^ • • • Ax^ exist as functions of 
/i, vanishing with A/i, and determined uniquely, so that the family is simple. 

Second, if the X/s are independent of t explicitly, the period also may be 
allowed to vary with /i, and this is quite conveniently accomplished by varying 
the initial instant t^. In this case if any determinant of the matrix 



(54) 



^x^ dx. 



ftc. 



dt 



(i,i = l.-.n), 



70 *'o 

does not vanish at ^ = ^^ + 27r, the increment of the initial constant correspond- 
ing to the omitted column may be considered arbitrary, and the other increments 
are then uniquely defined as functions of that increment and of A/i. For 
example, the determinant formed by omission of the nth column, on account of 
the linear relation (44), reduces to 



(55) 






70 



{t^l'*'»,j=l**'n — 1) 



'i 



where X^ is abbreviation f or X^ ( /a^ ; sCj^ , • • • , jb^ ) . Jtt does not enter explicitly 
in the X.% it is in general necessary to allow the period also to vary with fi, 
since in this case the determinant is always zero, as Poingab^ has pointed out.* 
This follows easily from the fact that the functions X.^ considered as functions 
of t along the initial periodic solution, where fAss fi^^ will themselves define a 

* MSth. Nauv., I, p. 90. 
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periodic solation of equations (42) ; a case which can happen for homogeneous 
linear equations with periodic coefficients only in the event that the equation 
in a>: 

(66) \e^(t, + 2ir)^8^w\^0 (»,i=l-n), 

has a root equal to unity, the notation being that of (28), (29), with ^^^ ss , so 
that in the present case 0^ is identical with dx./dxj^. 

Another illustration of Theorem IV of less obvious character refers to the 
equations considered in Theorem III. If there are k real positive numbers 
Xj, •• *, X^, such that 

(67) Lim JV^J»...JV^J*{X,-Xr} «0 

then similarly for the solutions 

(68) Lim JV^J' . . . J^i* { x^(t) ^x[{t)}=0 

the I's being in this case all equal to zero for i^ t^^ and consequently for all 
values of t. Other functions than powers of JVj, • • • , iV]^ might enter into sim- 
ilar relations. In particular, if the X's are all positive and greater than unity, 
the series defined by (80), (81), are absolutely convergent. 

§ 6. Integrals. Variation of parameters. 

We now consider the initial values ^^9 sb^^, • • • , sc^^ as parameters independent 
of the parameter r, and add to the preceding hypotheses the following, using the 
notation of § 8. 

VII. The functions X\ have continuous partial derivatives dX^jdXj^ which 
satisfy the conditions : 

if r is in R{€). 

If a similar condition were imposed with respect to every limit-point of R^ 
this would of course include I above, since the moduli of these derivatives would 
have an upper bound. 

Under these conditions the solutions x\ also have partial derivatives with 
respect to the initial parameters, satisfying the relations : 

(60) ^J + 2:x;(x.,.....«,,«/^-o 

parallel to (44). 

The existence of these relations has been interpreted by Bendixson and 
Lind£l5f f as the existence of n integrals a>p • - • , a>^ of the linear partial dif- 
ferential equation 

♦Floqubt, J. de I'iloole Norm., (2), 12, pp. 47-89. 
tLoo. oii 
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these integrals being determined by the condition that a>^ » Xj^ for ^ » ^j, where 
f J is an assigned value ; and the latter writer has specified in a simple way a 
realm of values (x^^ < ^^ x^^ t)^ within which those integrals may be known to 
exist. But for the present purpose it is necessary to inquire as to the depen- 
dence of that realm on the parameter r. 

If the value of t is considered fixed, the sole condition is that for the initial 
values %^ x^^y " 'f ^no^ ^^ solutions of equations (12') and (12'^) shall all be 
defined and lie within the realm D, over the closed interval t^' "t^j and since 
the moduli l^l^ |^^| have an upper bound independent of r it is possible to 
indicate such a realm independent of that parameter. The case ^^ > ^^ is easily 
treated in an analogous way. 

Let ^' be a point within the interval T^" ' T such that the smaller of the 
numbers t' — T^ and T^ if is 0^ and («{ , • . • , aj^) a point in D such that the 
lower bound of the sum S J ^^ — a;^ | is greater than the positive number S ; and 
let T be the smaller of the numbers 0^ SjiM; then the solutions are certainly 
defined and lie within D over the interval t^'-'t' provided | ^^ — ^ | = r and 
2 J sc^ ^ osj I ^ S/2 . With a change of notation which will be convenient later 
we have therefore the following theorem. 

Theobem V. If ^ is interior to the interval T^- •* T^ and (ajj , • • -, x'J) 
interior to the realm D, then l! is interior to an interval I within T^^ > T^ and 
(aj| , •••»*») i^^io^ to a realm A within D, such that for every t in /, every 
(Wj, • • •, x^) in A, and every r in H there exist the functions «j, . . ., o)^, 
Kf "*'^Nt (admitting continuous derivatives with respect to the arguments 
f, 0?,, • • •, se^, and satisfying the equations 

(61') ^+2:x,(«.,...,«.,o^-o, 

(61") -^ + 2:^;(«., •••.«., 0^-0, 

which are identities with respect to os^, • • *, a;^ m A, and t in /. The moduli 
of the partial derivatives have an upper hound independent of r, and the 
moduli of the Jacobians 

d{w^, '", a>J a(a)[, ..., <) 

d{x^, ..., 35^)' e(ajp '",xj' 

have a lower bound distinct from zero and independent of that parameter. 

As chosen above these integrals are uniquely defined by the conditions 
(o^ = a>][ = Xj^ for t s= t\ and are the solutions of equations (12'), (12''), at the 
instant t\ taking the initial values x^, • • *, a;^ at the instant t. 
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It seems needless to make a formal statement of a modified form of this 
theorem, analogous to Theorem I, whereby the existence of these integrals for 
equations (61'), for a certain realm A and a more extended interval for t than 
above, might imply the existence of the integrals of (61'') for the same realm 
and interval, but only for r in a certain subset Ii(€). The realm A and the 
set ^(e) could be readily specified by the use of the general inequality (17). 

The <o's have the same values for all points of any particular solution, and are 
of course the integrals or constants of integration, to be identified with the 
initial constants x^. To distinguish them from other sets of constants they may 
be called normal integrals or normal constants. The exact solution P(t) and 
the approximate solutions P''{t) are then defined respectively by the equations 

(62) a),.(ajp •••5«,5 = ««» ^'A^'n ""»K^^)'=^^^ 

where again for simplicity we suppose the initial constants independent of r. 
The results above show therefore that the equations 

a>. (aj^ + AflCj, . . ., x^ + Aaj^, t)^x^ + Aaj^, 

a>J(a;J + AajJ, •••,< + Aa?;, t)=^x^ + Ax„j 

can be solved so as to give the departures Ax.^ Ax[j of the solutions for any 
value f , in terms of the departures Ax^ of the initial constants, for solutions in 
the neighborhood of the assigned solution ; but it must be remembered that the 
restrictions on the magnitude of the increments Ax^ under which that solving 
can be known to be possible depend in general on the values of t and r, and it 
is conceivable that there should be no single neighborhood of the point 
(XjQ, "*f ^1,0) ^^^ which those restrictions are fulfilled simultaneously for all 
values of t and r. Inspection of the proof of the theorem on implicit functions 
above quoted shows, however, that those limitations on the increments are 
defined by certain inequalities to be satisfied by quotients of determinants, the 
denominators being the Jacobians just mentioned. Now for the interval of t 
and the range of initial values limited as specified in obtaining Theorem V, the 
moduli of the partial derivatives and of the reciprocals of the Jacobians have 
upper bounds independent of t and r. It is therefore possible to define the 
inequalities in question by the use of a set of positive numbers independent of t 
and r. Similar conclusions hold if the solutions P{t) and JP'(t) ^' cross" at 
any value in the interval named. 

For the following we require to consider the transformations of variables from 
(ajj, •••, a?^, t) to (©J, • ••, 0)^, t) or (©[, . .., ©j;, <), defined by the equations 

(68) a)^(a5i, '",x^,t)^w., o>J(«p --"i^n* «)**<• 

It is then shown by the foregoing that if (aSj^, • • •, xj) is interior to D and ^ to the 

interval T^-" T, then there is a closed set A^ of points (x^^+Ax^^ , • • • , x^+Ax^^) ^ 
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and a closed interval /^ or ^^ — • • • ^^ + 0, such that for any point on the solu- 
tion P(t) which starts in A^ for t = t^^ or any point on a solution I^{t) which 
crosses such a solution, the transformations are uniquely defined and have unique 
inverses. For example (a?^ . • •, x^) at t determines (a;^^^ • • -yX^) or (cd^, . . ., ©J 
which determines (»[ , • - • , «Jj) at t, which determines (asj^,, • • • , x^^) or (»[, • • •, ©;;)• 
There is a common realm within which there is a one-to-one correspondence 
between the sets of variables (a5j, •• ., sc^, ^), ((w^, • • •, o)^, t)^ (©[, .. ., w^, t). 
This realm contains in particular as interior points all points on the assigned 
solution P(^) for a certain ^interval. 

The method of variation of parameters may be defined as the employment of 
a transformation of variables such that the new variables are integrals or con- 
stants of integration for an approximate system of differential equations. If the 
process employed contemplates only a finite number of such transformations 
there is no question of convergence to be considered. But with respect to 
dynamical problems, like the problem of n bodies, this method, unless combined 
with others, may yield the required simplicity and symmetry of treatment at each 
stage only by involving an infinite sequence of transformations. In this case 
also it may be considered that one obtains an approximate solution by treating 
the variables at a certain stage of the process as constants, and then the question 
arises whether these approximate solutions converge toward the exact solution. 

The preceding theorems furnish means of testing this convergence for certain 
general classes of differential equations, for example where the approximating 
equations are defined by series as in § 4. For the determination of the exact 
solution for a certain interval by an infinite sequence of transformations, how- 
ever, one more condition must be satisfied ; the approximate solutions must he 
determine by the method for a certain common interval. As applied to trans- 
formations this would imply the condition that all the transformations employed 
must have a common realm of definition in the space where the coordinates are 
^1 * " ^n^ that realm having all points of the required solution as interior points. 
Under the hypotheses of the beginning of this section this condition is satisfied 
for the variation of parameters, if at each stage the integrals chosen are the 
normal integrals. 

Theorem VI. For differential equations of the type defined in theorem III 
there is a certain interval for t over which the solution starting at an interior 
point of D for t^^t^is the uniform limit of the approximate solutions furnished 
hy the method of variation of parameters^ if the integrals at each stage are 
chosen as the normal integrals. 

It is not essential that all approximate solutions have the same initial values 
x^ as the required solution, but simply such as lie in the neighborhood of that 
point and converge toward it. It is clear that many other choices of integrals 
might be made so as to satisfy the required conditions. 
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§ 7. Solomorphic functions. 

Attention has thus far been confined to real variables. It is known, however, 
that the method of difference-equations employed in the proofs applies equally 
well to analytic functions for complex values of the arguments, so that the results 
given above can be adapted to such cases also, with the interval ^o " ' ^o + ''^ 
replaced by a certain path of integration in the complex plane, leading to the 
point at which the value of the solution is sought. If the path of integration 
be defined by a real parameter it is su£Bicient to consider equations where the 
dependent variables, parameters, and functional values are complex, but the 
independent variable real. 

For equations of the type 

dz. 
(64) ^^^z,{z,,...,z^,ii,t) 

the linear equations of § 6 furnish a simple means of verifying the theorem of 
PoiNGAB^ relative to the analytic character of the solutions as functions of 
their initial values, or of the parameter ii. Let the functions Z^ be assumed 
complex, and holomorphic with respect to the complex arguments 2^9 • • '9 s^^^ ft) 
but simply continuous with regard to the real variable t . 

If the initial values s^io* " *9 ^^o ^^^ holomorphic with respect to fi^ihe solu- 
tions are. In particular the solutions are holomorphic with respect to the 
initial values.* 

For, let 



then equations (64) may be replaced by 

(66) -^ = X,{x, y, p, <r, t), "dt^^*^^' ^' f^' ^' ^)' 

which are equations in real variables, satisfying the conditions of § 6, if the 
initial values lie in the realm of regularity of the functions Z^. The derivatives 
dxjdp^ dxjda-y dyjdp^ ^yif^<^i therefore exist and are defined by the equations 

ax 

+ * 



dt dp ~ y^y dx^ dp df/j dp ) 

d dx^ ^ /dX.dxj dX^dyj\ 
di ~da^^\~dx^ 'da^'d^ ~da ) "*" 

dt dp~ ^\dxj dp '^ df/j dp J'^ 

dt d<T^y\dXj da '^ dy^ da ) "^ 



*Mkh, N<mv., I, p. 58. 
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By virtue of the conditions of monogeneity : 

^ ' dxj dy dy dx. dp dtr dtr dp 

We have therefore 

±(Sj^,djj\_^ldX,/dx,_dyj\dX,/dxjdyj\\ 

d(dx, dy\ I dX,(dx, dy,\ dXJdxj dy,\\ 

dt\8a'^dp )~'^\ dyj \dp dc )^ dx^ \diT ^ dp )]' 

a system of homogeneous linear equations satisfied by the 2n quantities 

^< _ ^i ^< , ^< 

dp 5<r ' da dp ' 

Because of the conditions on the initial values 

we have therefore for every value of t : 

A minimum value for the radius of convergence of the resulting power-series 
in p, can be given by reference to the general inequalities of § 8. 

§ 8. The canonical equations. 

We now consider briefly the special forms assumed by some of the preceding 
results for the canonical equations (1). The conditions imposed on the func- 
tions JC., X% will be satisfied if Fy F\ for all values of r in -R have all first 
and second derivatives with respect to as^ , * - * ) a;^ « y^ « " * 9 ^n ^ these derivatives 
being continuous functions of all arguments, including r in the sense that they 
satisfy the appropriate convergence conditions analogous to III of § 8, with 
respect to every limit-point of R. 

The linear equations of § 5 take the form* 

(71) ^^"rUy.^/^"^5y.5y/^; + lay,a,.)' 

dt ^\ dx,dx/j ^i^Vj'^O l^».-^/*r 



^PoiNCAB^ MHh. Nouv., I, p. 166. 
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which are also canonical, since they can be written 

^ ' dt~ dr,^' dt~~'di.- 

If we put 

(78) ?L*^^^'^''^^^>^''''"^*^^''*''^J 

and f , 7/ are derivatives with respect to /i. For derivatives with respect to the 
initial values the terms in { } are to be omitted. 

If the system in complex variables (64) is canonical, the corresponding real 
system is also canonical. For if we put 



a,, - x; + i/^ 1<, y, = y; + V- ly", F^ F'^- V- IF", 
equations (1) take the form 

^' ^ rf< " a^; ' dt — ai; ' dt " dx': ' dt "" By: ' 

by reason of the conditions of monogeneity. 

In employing transformations of variables it is convenient to choose such 
systems as retain the canonical form of the equations. It appears from the 
work of Jacobi* that if a new set of variables a^, • • •, a^, ^8^, • • •, ^8^, t^ i& 
defined by the solution of the equations 

(75) y.«a^. «. = 5^ 

where S{x^j • • • » a^^i i^n • • • » )8^ » ^^ *^ arbitrary function, then the trans- 
formed differential equations are 

For instance, if 8 satisfies the partial differential equation 

(77) -ar + ^'r"""'""' 5*;''"' ^'7='' 

and F is the som of F^ and F^, then equations (76) are 

da, dF, d^, dF, 



(78) 



dt 5/8, ' dt da. 



* Vorletungen Sber Dpumik, leotores SO, 36. 



